Information gain in tomography - A quantum signature of chaos 
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We find quantum signatures of classical chaos in various metrics of information gain in quantum 
tomography. We employ a quantum state estimator based on weak collective measurements of 
an ensemble of identically prepared systems. The tomographic measurement record consists of a 
sequence of expectation values of a Hermitian operator that evolves under repeated application of 
the Floquet map of the quantum kicked top. We find an increase in information gain and hence 
higher fidelities in the reconstruction algorithm when the chaoticity parameter map increases. The 
results are well predicted by random matrix theory. 
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Signatures of chaos in quantum mechanics appear in a 
variety of contexts. Examples include the level statistics 
associated with the random matrices of chaotic Hamil- 
tonians [1, 2], hypersensitivity of dynamics to pertur- 
bations [3, 4], the role of open quantum systems dy- 
namics (decoherence/measurement) in the emergence of 
chaos [5, 6], and the connection between chaos and the 
dynamical generation of entanglement [7-11]. In this let- 
ter we identify and analyze a new signature that unifies 
the characterization of chaos in both classical and quan- 
tum physics - chaos as a source for information gain in 
state estimation (tomography). 

At a fundamental level, chaos represents unpredictabil- 
ity, so this seems at odds with the goal of estimating 
an unknown state. On the flip side, however, this un- 
predictability represents the potential information to be 
gained in an estimation process. If everything is pre- 
dicted and known, we learn nothing new. The missing 
information in deterministic chaos is the initial condi- 
tion. In classical dynamics, a time history of a coarse 
grained trajectory at discrete times is an archive of infor- 
mation about the initial conditions given perfect knowl- 
edge about the dynamics. Moreover, if the dynamics is 
chaotic, the rate at which we learn information increases 
due to the rapid Lyapunov divergence of distinguishable 
trajectories. This information-theoretic picture is quan- 
tified by the Kolmogorov-Sinai (KS) entropy [12], which 
measures the information required to retrieve the increas- 
ingly fine-grained information about the initial condition 
that is necessary to maintain a constant coarse-grained 
prediction of the future chaotic trajectory. 

To probe the connection between quantum chaos and 
tomography, we consider a protocol based on the weak 
(nonprojective) collective measurement of an ensemble 
of Na identically prepared states that undergo well cho- 
sen dynamics [13, 14]. The time series of the measure- 
ment record provides the information used to reconstruct 
the initial condition. The dynamics is "informationally 



complete" if the time history contains information about 
an arbitrary initial condition. Such a protocol has been 
implemented through continuous time measurement of 
an atomic spin ensemble driven by external magnetic 
fields while it is monitored by a weakly coupled laser 
probe [15, 16]. This protocol allows us to explore how 
the chaotic nature of the dynamics is revealed in the in- 
formation content of the measurement record. 

Our goal in quantum tomography is to determine the 
state po (here, for a spin j) given an ensemble prepared 
in the state p® A ■ The system undergoes dynamics ac- 
cording to a prescribed unitary evolution, U(t), and is 
measured collectively as described above. The elements 
of the POVM, labeled by measurement outcomes X(t) at 
time t, are taken to be [17] 

E xw = ^=e Xp [-^(X(t)-J z (t)) 2 }- (1) 

We use the Heisenberg picture, so the observable being 
measured at time t is J z (t) = W{t)J z U(t), where J z is 
the z-projection of the collective angular momentum op- 
erator, J z — Y^iJz j an d ji^ is the projection for the i th 
spin. The Gaussian spread in the POVM elements, cr, is 
set by the shot noise of the probe. When the random- 
ness of the measurement outcomes is dominated by the 
quantum noise in the probe rather than the uncertainty 
AJ Z — \/Na Aj z ("projection noise"), quantum backac- 
tion is negligible, and the state remains approximately 
separable [17]. In this case the stochastic measurement 
record, normalized by the number of atoms, is well ap- 
proximated by 

M(i) = X(t)/N A = Tk(j'*(*)Po) + W{t), (2) 

where W{t) is Gaussian white noise with spread ct/Na- 
We parameterize the state in terms of a generalized 
Bloch vector r for a d-dimensional Hilbert space (d — 

d 2 — l ~ 

2j + 1) by p = I/d + Y, a =i r a-a, where form 
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an orthonormal basis of traceless Hermitian operators. 
Taking the record at discrete times M n = M(t n ), it then 
follows from Eq. (1) that in the weak backaction limit, 
the probability of finding measurement history M condi- 
tioned on the state r is 

p (M|r) <x exp ( " E ) ( 3 ) 
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where O na = Tr (j z {t n )E a ) and CT 1 = O t O is the 
inverse of the covariance matrix. The peak of this dis- 
tribution is the (unconstrained) maximum-likelihood es- 
timate of the Bloch vector, Yml = CO T M. The prob- 
lem of quantum tomography, therefore, reduces to linear 
stochastic state estimation. The eigenvalues of C _1 de- 
termine the relative signal-to-noise with which we have 
measured different observables (represented by its eigen- 
vectors). When the covariance matrix is full rank, the 
measurement record is "informationally complete." Due 
to finite signal-to-noise, however, the Bloch vector yml 
may not be associated with a physical density matrix 
with nonnegative eigenvalues, and therefore we must con- 
strain the solution. The final estimate f is found as the 
closest matrix according to 



r = argmin^[A/i - ^£> iQ r Q ] 2 , 



(4) 
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We find the solution efficiently via convex optimization. 

We now come to the central question of this work. How 
is the information content in the measurement record re- 
lated to the complexity of the dynamics encoded in U(t)7 
To simplify the analysis, we consider periodic application 
of a given Floquet operator U T , so that at the n th time 
step U(nr) — U™. The measurement record generated by 
such periodic evolution is generally not informationally 
complete; it lacks information about a matrix subspace 
of dimension > d— 2 out of the total dimension d 2 — 1 [18]. 
However, the condition that p be a positive matrix is a 
powerful constraint, that effectively allows "compressed 
sensing" [16, 19], resulting in high fidelity reconstruction 
with the available information in the measurement record 
generated by the orbit of a single U T . 

Our goal is to relate the information generating power 
of U T to the properties of quantum state reconstruction 
according to the protocol above. For this purpose, we 
consider a well-studied paradigm of quantum chaos: the 
kicked top [2] , described by the Floquet operator 



In our analysis, we fix a = 1.4 and choose A to be our 
chaoticity parameter. As we vary A from to 7, the clas- 
sical limit of the dynamics change from highly regular to 
completely chaotic. In the quantum description, as the 
dynamics becomes globally chaotic, and for j > 1, the 
Floquet operator has the properties of a random matrix 
picked from the appropriate ensemble. It is this ran- 
domness that leads to the analog of ergodic mixing for 
quantum systems. We take j = 10 in our studies, which 
is sufficient to achieve the statistics of random matrices, 
but small enough to be essentially quantum. 

We study the behavior of our reconstruction algorithm 
for an ensemble of 100 random pure states sampled from 
the Haar measure on SU(d), where d = 2j + 1 = 21. 
The dynamical evolution of the measurement record is 
generated by repeated application of the kicked top Flo- 
quet operator in Eq. (5). Figure la shows the fidelity 
of the state estimate, p, relative to the true state, l^o), 
•F = (V'oIpIV'o)j averaged over the ensemble of random 
states as a function of time, for different values of chaotic- 
ity in the kicked top. Two important features are ap- 
parent. As the level of chaos increases, both the rate 
of increase of fidelity and its value after 100 periods of 
control by the Floquet operator increases. When the 
classical description is globally chaotic, the correspond- 
ing quantum Floquet operator is well described as a 
random matrix and, together with the positivity con- 
straint, allows us to accomplish high fidelity quantum 
state reconstruction. We can further quantify the cor- 
relation between chaos and the performance of quan- 
tum state estimation using information theoretic met- 
rics. Brukner and Zeilinger [20] defined the informa- 
tion available in a quantum measurement of a system, 
E, by the uncertainty of the outcomes summed over a 
set of mutually complementary experiments. Rehacek 
and Hradil showed that this uncertainty is equal to the 
Hilbert-Schmidt distance between the true and estimated 
state in quantum tomography, averaged over many runs 
of the estimator, E = (Tr{(po — p) 2 }) [21], which can 
be expressed as the total uncertainty in the Bloch vec- 
tor components, E = X)a((^ r «) 2 ) • These uncertain- 
ties are always greater than the Cramer-Rao bound, 
((Ar Q ) 2 ) > [F" 1 ]^, where F is the Fisher information 
matrix associated with the conditional probability distri- 
bution, Eq. (3), and thus E > Tr F _1 . In the limit of 
negligible quantum backaction, we saturate this bound, 
as our probability distribution is Gaussian, regardless of 
the state. In that case, the Fisher information matrix 
equals the inverse of the covariance matrix, F = C" 1 , in 
units of N\/a 2 . Thus, a metric for the total information 
gained in tomography is the inverse of this uncertainty, 
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(5) which measures the total Fisher information. Quantum 
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FIG. 1. Information gain in tomography as quantified by various metrics, given a measurement record generated by iterations 
of a quantum-kicked-top Floquet map U T = exp{— iAjf/(2j)}exp{— iaj x }, for a spin j = 10. The value of a is fixed and A 
serves as the chaoticity parameter, varying from regular dynamics, A = 0.5, to fully chaotic, A = 7.0. All results are shown as 
the average of 100 Haar-random pure states, (a) Fidelity of state reconstruction, (b) The Fisher information of estimating the 
parameters that define the state, (c) The Shannon entropy of the normalized eigenvalues of the inverse of covariance matrix 
of the likelihood function. The maximum possible Shannon entropy is shown as the dashed-dotted line. In all cases, both the 
rate of growth and the final value of the information metric are increased with higher values of the chaoticity parameter, A. 
In the fully chaotic regime, A = 7.0, the results are well predicted by tomography performed with a the measurement record 
generated by a Haar-random matrix picked from the circular orthogonal ensemble (COE). The COE results, averaged over 100 
random pure states, are plotted as the dashed line in (a)-(c). 



tomography can also be viewed as a form of "parameter 
estimation," i.e., estimate the Bloch vector components 
that define po- The Fisher information, then, quantifies 
how well our estimator can predict these parameters from 
the data, regardless of the state. 

In Fig. lb, we plot J as a function of time, gener- 
ated by repeated application of the kicked top dynamics 
described above. As before, we see the close correlation 
between the degree of chaoticity and the information gain 
in tomography. Note that the inverse covariance matrix 
is never full rank in this protocol, and J is always ill 
defined. We rectify this, regularizing C _1 by adding to 
it a small fraction of the identity matrix, similar to the 
Tikhonov regularization (see, e.g., [22]). In this way, our 
estimator ignores the directions in the space of observ- 
ables that are largely unmeasured, and then makes the 
best guess consistent with the positivity constraint. 

There is a close relationship between Fisher Informa- 
tion and fidelity as a metric for information gain. When 
Po is a pure state, the average Hilbert-Schmidt distance, 
E = 1/J = l-(Trp 2 )-2(J") [23]. A correlation between 
chaos in the dynamics and the information gain as seen 
in the average fidelity implies that the Fisher information 
will exhibit the same correlation. Moreover, the Fisher 
information can be further related to a true information 
metric - the mutual information I[r; M] - defined as the 
information we obtain about r from measurement record 
M, which is given by J[r; M] = H(M) - H(M\r) [24]. 
Here, H is the Shannon entropy of the given probabil- 
ity distribution. Assuming perfect knowledge of the dy- 
namics, the entropy of the measurement record, i?(M), 
irrespective of the state, is due solely to shot noise, and 



is thus constant. This is analogous, in the classical case, 
to the entropy associated with equal a priori probabil- 
ity to find a trajectory in one of the coarse grains of 
phase space. Neglecting irrelevant constants, the mutual 
information between the Bloch vector and a given mea- 
surement record is then specified by the entropy of the 
conditional probability distribution, Eq. (3), 

X[r; M] = -H(M\r) = -~ log (detC) = log(l/y), (7) 

where V is the volume of the error-ellipsoid whose semi- 
major axes are defined by the covariance matrix. 

In order to maximize the information gain, we seek 
the dynamics that maximizes 1/V = ^/det (C _1 ). An 
important constraint is that after time t n , 

Tr(C- 1 ) = ^(a, Q ) 2 = n||O(0)|| 2 , (8) 

i,a. 

where ||O(0)|| 2 = £ Q Tr (0(O)S Q ) 2 is the Euclidean 
square norm, with O(0) = j z for our example. The 
right-hand-side of this equation is independent of U and 
increases monotonically with time. It then follows from 
the inequality of arithmetic and geometric means, 

det^- 1 )^ ^Tr(C- 1 )^=(^||O(0)|| 2 ) D . (9) 

where D = d 2 — 1 is the rank of the regularized covariance 
matrix. The maximum possible value of the mutual in- 
formation is attained when all eigenvalues are equal and 
the above inequality is saturated, implying that the error 
ellipsoid is a hypersphere. At a given time step, the dy- 
namics that gives the largest mutual information is the 
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one that mixes the eigenvalues most evenly. We quan- 
tify this by the Shannon entropy H of the eigenvalues of 
the inverse of covariance matrix, A Q , normalized to as a 
probability distribution, p a = A Q /Tr (C _1 ). Figure lc, 
shows H{p) as a function of time generated by the kicked 
top Floquet map. Again, we see a close relationship be- 
tween the degree of chaoticity of the map and this metric 
for information gain in tomography. This makes sense 
physically. In order to extract the maximum information 
about a random state, we must measure all components 
of the Bloch vector with maximum precision. Given fi- 
nite time, we obtain the best estimate by dividing equally 
between all observables. 

Finally, we show that the information gain generated 
by the quantum-chaotic dynamics is fully consistent with 
random matrix theory. The standard kicked top dynam- 
ics is time-reversal invariant without Kramer's degener- 
acy [2]. For parameters in which the classical dynamics 
is globally chaotic, we thus expect the Floquet opera- 
tor to have the statistical properties of a random matrix 
chosen from the circular orthogonal ensemble (COE) [2]. 
Figure 1 shows the behavior of the fidelity, the Fisher 
information, and Shannon entropy of the inverse of the 
covariance matrix as a function of the number of appli- 
cations for a measurement record generated by a typi- 
cal random unitary picked from the COE (dotted line). 
We see excellent agreement between our predictions from 
random matrix theory and the calculation for the evolu- 
tion by the kicked-top Floquet map in completely chaotic 
regime, a = 1.4, A = 7. Note that the maximum Shannon 
entropy, H max = log(d 2 — 1), is never attained. This is 
because our measurement record is generated by an orbit 
of repeated application of a single U T . Nonetheless, the 
measurement record, together with the positivity con- 
straint, is sufficient to yield high fidelity reconstruction 
when the Floquet map is deep in the chaotic regime [18]. 

In summary, complex dynamics reveals more informa- 
tion about the initial condition as one observes the sys- 
tem over the course of time. Classically, chaotic dynam- 
ics, with its Lyapunov sensitivity to the initial condi- 
tions, generates an exponentially expanding archive of 
information about the initial state if we are able to track 
a trajectory with a constant coarse-grained resolution. 
Similarly, we found that the rate at which one obtains 
information about an initially unknown quantum state 
in quantum tomography is correlated with the extent of 
"quantum chaos" in the system, i.e., the degree to which 
the unitary dynamics maps a localized coherent state to 
a random state in Hilbert space. Quantum tomography 
thus provides a forum in which to unify the the notions of 
complexity of chaotic dynamics in classical and quantum 
worlds. We quantified the information gain in a variety of 
metrics, including reconstruction fidelity, Fisher informa- 
tion, mutual information, and Shannon entropy. When 
the system is fully chaotic, the rate of information gain 



is well predicted by random matrix theory. This novel 
signature of chaos can be explored using current experi- 
mental techniques in the setting of cold atoms interacting 
with lasers and magnetic fields. 
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